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SCALE-INVARIANT  THEORY  OF 
OPTICAL  PROPERTIES  OF  FRACTAL  CLUSTERS 

VADIM  A.  MARKEL'.  LEONID  S.  MURATOV. 

MARK  I.  STOCKMAN z:  ana  THOMAS  F  GEORGE : 

Institute  of  Automation  &  Electrometry.  Siberian  Branch  of  the  USSR  Academy  of 
Sciences 

•  Departments  of  Physics  6c  Astronomy  and  Chemistry.  Center  for  Electronic  6:  Electro- 
optic  Materials.  239  Fronczak  Hall.  State  University  of  New  York  at  Buffalo.  Buffalo. 
NY  14260 

INTRODUCTION 

'"'r  Fractals  are  scale  self-similar  mathematical  objects  possessing  nontrivial  geometri¬ 
cal  properties^!  1.21*.  There  exist  various  physical  realizations  of  fractals  3-5..  and  here 
we  shall  consider  optical  properties  of  what  we  believe  to  be  one  of  the  most  important 
suiywralizations.  namely,  fractal  ciusters^to  be  named  below  simpiv  as  fractals).  . 

A  fractal  cluster  is  a  system  of  interacting  material  particles  called  monomers.  This 
system  is  self-similar  (in  a  statistical  sense)  with  respect  to  scale  transformations  in  an 
intermediate  region  of  sizes  r.  Rq  <c  r  <  Rc>  where  Ro  is  a  characteristic  separation 
between  nearest  monomers,  and  Rt  is  the  cluster  total  size. 

A  consequence  of  the  self-similarity  is  the  scaling  (power-like)  form  of  the  pair 
(density-density)  correlation  function  g(r),  mean  number  N  of  the  monomers  in  the  the 
fractal  and  its  mean  density  p, 

g(r)  =  (D/4nR*)(r/R0)D-3  .  N  -  (r/Rof  .  p  ~  {Rc/R0)D~i  ■  (1) 

where  the  index  D  is  the  fractal  (external  Haussdorff)  dimension.  A  fractal  is  called 
nontrivial  if  D  <  3.  As  one  can  see  from  Eq.  (1),  the  mean  density  of  a  nontrivial  fractal 
is  asymptotically  (for  R<.  2>  Ro)  zero.  This  feature,  together  with  strong  pair  correlation, 
is  the  reason  for  the  great  role  that  the  density  fluctuations  play  in  the  fractals.  Below 
we  shall  show  that  large  correlated  fluctuations  of  local  fields  play  a  decisive  role  in 
determining  optical  properties  of  fractals.  In  particular,  there  exists  scaling  of  the  optical 
absorption  and  the  density  of  the  dipole  excitations  (eigenstates)  of  the  fractal. 

The  linear  optical  (dipole)  polarizabilities  of  fractals  have  been  studied  theoretically 
6-91.  However,  the  scale-invariant  theory,  which  is  the  subject  of  the  present  report,  has 
not  been  published  in  the  literature. 

FRACTAL  OPTICAL  PROPERTIES  AND  SCALING 

Let  us  consider  a  model  of  a  cluster  as  a  system  consisting  of  .V  polarizable  particle 
'called  monomers)  located  at  the  points  r,,  i  =  1....V.  with  the  dipole  interaction 
between  them,  and  Rc  <  A.  where  A  is  the  light  wavelength.  The  light-induced  dipole 
(transition)  moments  d,  on  monomers  obey  the  well-known  system  of  equations 

d,a  =  XoEa0)  ~  Xo  V  '  \8as  -  3 n^;)ny-"  1  r~%d.3  .  12) 

:  =  1 

where  the  Greek  subscripts  stand  for  tensor  components  ( summation  over  repeated  in¬ 
dices  is  implied),  and  the  Latin  indices  stand  for  ordinal  numbers  of  monomers:  xo  IS  the 


dipolar  polarizability  of  an  individual  (isolated!  monomer:  E'01  is  the  electric  field  of  the 
external  (exciting;  radiation:  ri;  =  r.  -  r,:  and  nUJI  —  r,,,  r... 

The  fractal  specificity  of  the  considered  problem  is  due  to  the  divergence  at  small  r 
in  the  mean- field  analog  of  Eq.(2).  E(0>  -  \0  J  d  -  3r(rd)/>:  r-3^(r ;dr  (cf.  Refs.  ; 8.9] ). 
This  divergence  is  an  indication  of  the  dominating  role  of  large  fluctuations  and  existence 
of  scaling  in  fractals.  In  this  brief  report  we  can  only  outline  the  theory  and  then  briefly 
discuss  the  results  obtained. 

The  main  object  of  the  theory  is  the  optical  (dipole)  polarizability  \(A’)  of  the 
fractal  (per  one  monomer)  as  a  function  of  the  variable  X  =  Rey^1.  The  polarizability 
is  formally  found  as  a  decomposition  in  terms  of  the  exact  eigenvectors  and  eigenvalues  of 
Eq.  (2).  L  sing  this  decomposition,  the  exact  sum  rules  for  the  polarizability  are  obtained 


-  /  lmX„3(X)dX  =  6ad  .  P  /  ReX„a(X)dX  =  0  ,  P  /  XImy(X)dX  =  0  . 

(3) 

where  P  denotes  the  principal  vaiue  of  integral. 

The  next  step  is  to  assume  the  scaling  of  dipole  excitations  of  the  fractal.  Then  the 
•.unction  Imy(.Y)  can  be  shown  to  have  the  scaling  form  (note  that  Re\  does  not  scale), 
if  the  resonance  in  the  monomer  has  a  high  Q-factor.  Q  =  f^Imy,)’1  -1.  From  (3)  it 
follows  that  Imy  scales  as  the  function  of  iJYl,  thus  having  the  form 

.  (4) 

where  the  index  d0  is  called  the  optical  spectral  dimension,  0  <  d0  <  1.  Then,  using  the 
fact  that  Imy  describes  absorption  and  the  total  absorption  of  a  fractal  in  the  scaling 
region  should  not  depend  on  choice  of  the  minimum  length  Ro,  from  Eq.  (4)  we  obtain 
the  scaling  transformation  law  iX|  oc  1  \  Using  this,  one  can  make  sure 

that  it  possible  to  construct  only  one  quantity  Lx  which  has  dimensionality  of  length 
and  is  invariant  under  the  scaling  transformation, 


Ro(Rl\X\) 


(do  —  1  )/(3  —  D) 


This  quantity,  in  accord  with  Alexanders  strong  localization  hypothesis  [10],  should 
coincide  with  the  coherence  length  of  the  fractal  excitations.  Equation  (5)  defines,  thus, 
the  dispersion  law  for  the  dipole  excitations  of  the  fractal,  with  1/Lx  playing  the  role  of 
wavevector.  The  conditions  of  scaling  Ro  Lx  <  Rc  and  X  with  the  aid 

of  Eqs.  (5)  and  (1)  axe  rewritten  as 

R;3Q-\  v^/o-ri/U-^)  <  x  <  R//3  .  (6) 

Hence,  if  the  necessary  condition  Q  >  1  is  met.  scaling  occurs  at  small  I  AT.  In  the  vicinity 
of  the  frequency  u>o  of  the  resonance  in  the  monomer,  the  variable  A*  is  proportional  to 
detuning  Cl  =  u>  -u>0,  where  u>  is  the  light  frequency.  Hence,  the  scaling  region  does  exist, 
occupying  the  center  of  the  fractal  absorption  band.  Equation  (6)  can  be  considered  as 
a  self-consistent  validation  of  the  scaling  assumption. 

Using  the  above  mentioned  exact  decomposition  of  \  and  Eqs.  (3),  one  can  show 
for  the  scaling  region  (6)  that  the  following  relation  exists  between  Imy  and  the  density 
v{X)  of  the  dipole  excitations  (eigenstates)  of  the  fractal: 

Irnx(X)  =  =<  .  (7) 


Thus,  the  density  of  states  and  optical  absorption  scaie  with  the  same  exponent  d0  —  1 .  and 
the  index  do  is  the  counterpart  for  the  ciipoie  excitations  of  the  vibration  i  fracton )  spectral 

dimensionality  d  introduced  by  Alexander  and  Orbach  11..  The  dipole  excitations  tend 
to  surface  plasmons  in  the  trivial  (D  —  31  limit,  while  vibrations  tend  to  phonons.  The 
surface  plasmons  possess  a  spectral  gap.  while  phonons  do  not.  which  is  the  reason  why 
the  dispersion  law  i  5)  differs  from  us  counterpart  10.11  .  For  the  same  reason,  the  index 

d0  does  not  numerically  coincide  with  d  |  cf.  0  <  da  <  1.  while  0  <  d  <  3). 

The  failure  of  the  mean-field  approach  and  existence  of  scaling  imply  that  fluctu¬ 
ations  of  the  local  fields  in  the  fractals  play  a  decisive  roie  in  determining  their  optical 
properties.  From  the  above  mentioned  expansion  of  \  in  the  terms  of  the  fractal  eigen¬ 
states.  it  is  possible  to  obtain  the  exact  relation  (an  analog  of  the  optical  theorem) 
describing  the  relative  dispersion  G  of  the  local  field  E,  acting  upon  the  ith  monomer  in 
the  fractal. 

G  =  (  E«i2)  .  E(0)|2  =  Q~l  (1  -  Q2X2)  Im\.s  .  (8) 

In  the  scaling  region  (6).  from  Eq.  (8)  it  follows  that 

G~~QRlXd°~\  (9) 

Thus,  the  fluctuations  are  large:  the  dispersion  of  E,  is  proportional  to  Q  »  I.  reaching 
its  maximum  at  the  wings  of  the  scaling  region,  i.e..  for  .  A'  <  1. 

NUMERICAL  SIMULATION. 

Three  types  of  fractals  have  been  generated  by  the  Monte-Carlo  method:  random 
walks  (D  =  2),  self-avoiding  random  walks  ( D  =x  1.7)  and  the  Witten-Sander  clusters 
121  ( D  %  2.5).  The  mean  number  of  monomers  in  the  cluster  was  i V  =  25  —  50,  and 
the  number  of  clusters  was  from  several  hundred  to  thousand.  For  each  cluster  the 
fundamental  system  (2)  was  solved,  and  averaging  over  different  fractals  was  performed. 
The  numerical  results  show  that  the  absorption  and  density  of  states  obey  the  seeding 
prediction  (4)  and  (7).  The  example  of  such  results  is  presented  in  Fig.  1.  The  the  values 
found  for  the  d0-exponent  from  the  figure  and  similar  data  for  other  fractals  are  given  in 
Table  1. 


rig.  i 

Double-logarithmic 
plot  of  lm\(A”)  calcu¬ 
lated  directly  (circles) 
and  expressed  through 
the  eigenstate  density 
with  the  aid  of  Eq.  (7) 
(squares').  The  stra¬ 
ight  lines  are  found  by 
the  linear  regression. 


Table  1 


Fractal  type:  Self-avoiding  walks  Random  walks  Witten-Sander's 
d0  found  from  Im\:  0.43  =  0.03  0.38  =  0.03  0.49  -  0.02 

do  from  eigenstate  density:  0.54  —  0.08  0.33  —  0.05  0.51  —  0.06 

CONCLUDING  DISCUSSION 

The  scaling  of  optical  absorption  in  nontrivial  fractals  is  suggested  bv  the  failure  of 
the  mean-field  approach.  Assuming  scaling  (4).  we  have  finally  found  the  condition  (6) 
of  it.  which  means  that  the  scaling  assumption  is  self-consistent.  Scaling  is  confirmed  by 
numerical  simulation,  and  the  values  of  d0  are  found  for  three  types  of  fractals.  Also, 
scaling  is  intimately  coupled  to  large  fluctuation  of  local  fields  in  nontrivial  fractals.  Such 
fluctuations,  as  shown  in  the  framework  of  the  binary  approximation  13!.  bring  about 
giant  enhancement  of  the  nonlinear  parametric  scattering  of  light  from  fractals.  The 
selective  photomodification  of  fractals  is  also  due  to  these  fluctuations  14; .  The  features 
mentioned  indicate  the  fractals  to  be  promising  media  for  nonlinear  optics  and  optical 
information  processing. 
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